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ABS1BAQX 


The  poor  aerodynamic  design  of  the  superstructures  of 
today's  Navy  ships  create  a  highly  turbulent  airwake  that  can 
make  shipboard  helicopter  operations  quite  hazardous.  This 
study  is  part  of  a  longer-term  project  to  tailor  airflow  over 
a  ship's  helicopter  deck  in  order  to  improve  the  poor  quality 
of  the  airflow.  This  airflow  is  thought  to  be  largely 
responsible  for  several  costly  blade  strike  mishaps  involving 
H-4  6  helicopters  during  start-up  and  shut-down.  Numerous 
computer  simulations  were  conducted  using  the  "Phoenics" 
Computational  Fluid  Dynamics  code  to  simulate  airflow  over  a 
backward-facing-step.  The  latter  represents  a  simple  2-d 
model  of  flow  behind  a  typical  hangar/flight  deck  combina¬ 
tion.  In  each  run  a  deflector  of  different  size,  orientation 
or  porosity  was  placed  in  a  specific  location  in  an  attempt  to 
reduce  the  size  of  the  recirculation  zone,  the  velocities  and 
the  turbulence  levels.  Of  the  studies  involved,  a  vertical 
deflector  offset  downstream  with  its  lower  edge  at  the  top  of 
the  step  produced  the  best  overall  results. 
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I.  INTRODUCTION 


The  Navy's  ability  to  conduct  helicopter  operations  is 
more  important  today  than  ever  before.  Helicopters  are  a 
major  contributor  to  the  Antisubmarine  and  Antisurface  Warfare 
mission  areas  and  are  equally  important  in  the  role  of 
material  supply  between  fleet  ships,  both  underway  and  at 
anchor.  The  twin  rotor  helicopter,  the  H-46,  is  commonly  used 
in  this  resupply  mission  because  of  its  large  load  carrying 
ability  and  tolerance  to  varying  wind  conditions  while 
hovering  over  ships. 

The  ability  of  the  H-46  to  operate  in  various  weather  and 
wind  conditions  is  vital  to  the  smooth  and  efficient  operation 
of  the  fleet.  However,  due  to  an  increasing  number  of  rotor 
strike  mishaps,  the  overall  availability  of  the  H-46  has  been 
reduced.  Most  of  these  blade  strike  mishaps  have  occurred  on 
AOR,  LPA  AND  LPH  type  vessels  (Figure  1)  ,  during  rotor 
engagement  or  rotor  shutdown,  with  winds  of  at  least  moderate 
strength. 

It  is  generally  thought  that  the  main  reason  for  most 
blade  strike  mishaps  is  the  turbulent  wind  conditions 
encountered  on  the  ship's  helicopter  deck.  The  only  presently 
known  solution  to  this  problem  is  a  reduction  of  the  rotor 
engage/disengage  envelope.  The  latter  prescribes  allowable 
safe  operation  limits  based  on  the  wind  direction  and  speed. 
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Figure  1.  AOR  Class  Ship  [Ref.  1] 


The  highly  turbulent  conditions  often  found  on  the 
helicopter  deck  are  a  result  of  the  interaction  between  the 
wind  and  the  ship's  superstructure.  When  Navy  ships,  were 
first  designed,  their  mission  dependence  on  helicopters  was 
not  foreseen  and  so  little  attention  was  paid  to  the 
aerodynamics  of  the  superstructure.  Unfortunately,  there  has 
been  little,  if  any,  improvement  in  superstructure  design, 
with  regards  to  airflow,  on  today's  newer  ships.  Thus,  the 
environment  around  the  ships,  in  which  the  helicopter  must 
operate,  remains  quite  hazardous. 

The  Naval  Postgraduate  School  (NPS)  has  undertaken  a 
series  of  studies  on  ship  airwake  tailoring,  with  the  goal  of 
reducing  the  level  of  turbulence  and  the  size  of  the  separated 
zones  on  the  helicopter  deck.  This  reduction  will  result  in 
a  larger  engage/disengage  envelope,  which  in  turn,  will 
increase  the  percentage  of  time  that  the  H-46  can  safely 
operate.  This  particular  study  is  the  first  in  a  series  at 
NPS  to  attempt  to  achieve  this  goal  and  will  consider  only  the 
feasibility  of  such  endeavors  by  tailoring  a  two-dimensional 
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flow  over  a  backward  facing  step.  The  "Phoenics"  computation¬ 
al  fluid  dynamics  ( CFD)  code  will  be  utilized  to  simulate  the 
two-dimensional  flow  and  the  tailoring  of  the  airflow  will  be 
accomplished  in  the  code  with  strategically  placed  deflectors. 

First,  a  more  in-depth  look  at  the  role  of  the  H-46  in 
fleet  operations  will  be  addressed;  this  will  be  followed  by 
a  short  review  of  bluff  body  aerodynamics,  ship  airwakes,  and 
any  wake  modification  methods  used  up  to  the  present.  After 
the  computer  simulations  are  made,  they  will  be  analyzed  and 
recommendations  made  on  step/deflector  configurations  for  the 
next  part  of  th^  program,  which  will  be  an  attempt  at 
experimental  verification  of  the  tailoring  of  two-dimensional 
flow  over  a  backward  facing  step  in  the  NPS  low  speed  wind 
tunnel . 
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II.  BACKGROUND 


A.  HELICOPTER  OPERATIONS 

Navy  fleets,  which  must  travel  worldwide,  cannot  hope  to 
accomplish  their  mission  without  the  "Mobile  Logistic  Support 
Force"  (MLSF)  ships.  [Ref.  2:p.  1.1]  This  group  of  ships 
consists  of  the  following  ship  classes:  AE,  AO,  AOR,  AFS ,  and 
AOE .  These  ships  sail  with  the  fleet  carrying  most  of  the 
supplies  required  for  that  particular  deployment,  such  as 
food,  fuel,  and  material  goods.  There  are  two  basic  methods 
to  transfer  these  supplies  from  ship  to  ship,  connected 
replenishment  (CONREP)  and  vertical  replenishment  (VERTREP) . 
[Ref.  2 : p.  1.3] 

CONREP  is  a  method  that  involves  two  ships  steaming  side 
by  side  within  80  to  200  feet  of  each  other  and  transferring 
supplies  by  means  of  cables  strung  between  them.  This  method 
is  commonly  used  for  refueling  ships  and  transferring  loads 
which  are  too  heavy  for  VERTREP. 

The  preferred  method  of  replenishment  is  with  helicopters, 
a  method  more  commonly  referred  to  as  VERTREP.  Supply  items 
are  placed  in  large  cargo  nets,  lifted  by  the  helicopters,  and 
transferred  to  the  appropriate  ships.  The  specific  advantages 
of  VERTREP  are: 

(1)  Reduction  in  time  required  to  replenish  the  supported 
forces  or  units. 
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(2)  Reduction  or  elimination  of  time  that  screening  ships 
are  off  station. 

(3)  Reduction  of  the  number  of  personnel  involved. 

(4)  Capability  of  replenishing  units  in  a  dispersed 
formation. 

(5)  Capability  of  replenishing  units  engaged  in  tasks  which 
make  it  impossible  for  them  to  come  along  side. 

(6)  Capability  of  replenishing  units  in  heavy  weather 
conditions  when  alongside  steaming  is  hazardous  or 
impossible. 

(7)  Capability  of  replenishing  units  on  station  in  shallow 
water  or  at  anchor. 

As  stated  before,  the  twin  rotor  H-4  6  has  become  the 
helicopter  of  choice  for  this  mission  and  its  continuous 
availability  in  various  weather  and  wind  conditions  has  become 
vital  in  fleet  operations.  Blade  strikes,  which  are  mishaps 
that  occur  when  a  turning  rotor  blade  impacts  the  helicopter 
fuselage,  cause  damage  that  must  be  repaired  before  the 
helicopter  is  available  again  and  are  numerous  enough  to 
threaten  a  reduction  in  the  safe  operating  envelope  of  the  H- 
46. 

Though  blade  strikes  typically  occur  at  a  very  low  rotor 
RPM  (usually  at  the  start  of  rotor  engagement  or  the  end  of 
rotor  shutdown,  when  rotor  RPM  is  about  20%  of  normal),  they 
still  pose  a  significant  danger  to  the  aircrews  and  ground 
personnel  as  well  as  the  aircraft.  Through  the  middle  of  1989 
there  have  been  at  least  75  such  incidents  aboard  ships 
ranging  from  little  or  no  damage  to  complete  loss  of  the 
airframe.  The  dollar  cost  can  range  from  just  man  hour  costs 
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when  only  inspections  are  required  for  a  minor  strike,  to 
upwards  of  $500,000  for  a  blade  strike  which  involves  a  sudden 
stoppage  to  the  drivetrain  system.  [Ref.  3:p.  4]  If  the 
entire  airframe  is  lost,  the  cost  can  escalate  to  well  over 
$2.7  million,  which  is  the  approximate  cost  of  a  replacement 
H-4  6 . 

Hidden  costs  are  somewhat  harder  to  measure,  but  are 
equally  costly.  An  interruption  in  the  UNREP  operation  can 
lead  to  changes  in  the  entire  fleet  schedule.  This  may  cause 
numerous  other  operations  to  be  changed  or  canceled  to  allow 
the  resupply  mission  to  continue. 

To  avoid  blade  strikes  while  operating  on  shore,  or  from 
a  ship,  the  helicopter  crews  use  a  chart  which  is  essentially 
a  go/no  go  chart  which  tells  them  if  wind  direction  and  speed 
will  allow  a  safe  rotor  engagement  or  disengagement.  The 
generic  envelope  for  the  H-46  is  shown  in  Figure  2.  There  are 
also  ship  specific  envelopes,  an  example  of  which  is  shown  in 
Figure  3. 

These  envelopes  are  developed  through  dynamic  interface 
testing  done  by  the  Naval  Air  Test  Center  (NATC)  .  It  is  a 
long,  laborious,  and  expensive  process  which  is  valid  only  for 
the  particular  combination  of  ship/helicopter  being  tested. 
A  description  of  the  testing  process  is  given  by  Madey  and 
Whitmer  [Ref.  5].  Unfortunately,  due  to  the  extraordinary 
variability  of  wind  and  sea  state  conditions,  it  is  nearly 
impossible  to  document  a  completely  safe  operating  envelope. 
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Figure  2.  Generic  Operating  Envelope  [Ref.  4] 


This  is  certainly  true  for  the  H-46  and  the  AOR  ship 
combination  in  which  a  significant  number  of  the  blade  strike 
incidents  have  occurred  while  operating  inside  of  the  rotor 
engage/disengage  envelope. 


Figure  3 .  Daytime  Operating  Envelope  for 
BB-63  [Ref.  4] 


This  predicament  has  led  to  the  suggestion  that  the 
problem  may  be  solved  by  simulation  [Ref.  6:p.  2).  To  achieve 
this  simulation  accurately  there  is  a  need  to  predict  the 
freestream  airflow  over  the  ship,  the  ship  motion,  and  the 
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motion  of  the  helicopter.  Healey  [Ref.  6:pp.  14-58]  looks  at 
what  has  been  done  in  these  fields  and  lists  several  other 
references.  The  program  underway  at  NPS  is  attempting  to  make 
detailed  air-wake  maps  of  model  ships  for  scaling  to  full- 
size.  So  far,  visualization  of  the  flow  around  a  model  of  a 
DD-963  class  destroyer  has  been  completed  [Ref.  7]  and  a 
similar  study  of  an  AOR  class  ship  is  presently  in  progress. 

Until  that  time  when  a  completely  accurate  rotor  engage/ 
disengage  envelope  is  developed,  other  alternatives  must  be 
looked  at.  This  paper  explores  one  such  alternative: 
tailoring  the  ship  airwake  in  order  to  reduce  the  size  of  the 
separation  zones  and  the  levels  of  turbulence  encountered  on 
the  helicopter  deck.  This  approach  will  not  only  eventually 
lead  to  a  safer  operating  area  for  the  helicopters  but,  if 
successful,  will  increase  the  size  of  the  rotor  engage/ 
disengage  envelope,  thus  allowing  the  helicopters  to  operate 
a  greater  percentage  of  time. 

B.  BLUFF  BODY  AERODYNAMICS 

Though  the  experimental  part  of  this  paper  will  be 
conducted  as  a  two-dimensional  problem,  it  is  an  important 
first  step  to  understanding  the  three-dimensional  airflow 
experienced  by  the  ship,  the  helicopter  deck,  and  subsequently 
the  helicopter.  Bluff  body  aerodynamics  and  their  relation¬ 
ship  to  this  problem  will  first  be  examined.  This  is  not 
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meant  to  be  a  detailed  discourse  on  the  subject,  as  the 
references  to  be  cited  more  than  adequately  cover  the  subject. 

By  definition,  a  bluff  body  is  one  in  which,  for  given 
flow  conditions,  there  is  a  massive  separated  region  in  its 
wake.  It  becomes  apparent,  after  first  observing  the  design 
of  any  large  class  Navy  ship  (Figure  1)  and  then  a  schematic 
of  the  observed  flow  over  a  model  ship  (Figure  4) ,  that  Navy 
ships  indeed  can  be  considered  three-dimensional  bluff  bodies. 


To  start  to  appreciate  the  complexity  of  the  flow  around 
a  bluff  body,  one  only  has  to  look  at  a  study  by  Hunt,  Abell 
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Peterka,  and  Woo.  [Ref.  8:pp.  179-200]  They  detected  the 
presence  of  an  inverted  U-shaped  vortex,  whose  ends  remained 
in  contact  with  the  ground,  on  the  downwind  side  of  the  body 
and  numerous  horseshoe  vortices  that  wrap  themselves  about  the 
upstream  base  of  the  body  and  trail  downstream  (Figure  5) .  In 
addition,  turbulence  causes  the  reattachment  point  to  be 
highly  unstable,  and  alters  the  flow  field  around  the  body  by 
producing  increased  mixing  near  the  separated  shear  layers 
(Ref.  9]. 


Figure  5.  Mean  Streamline  Patterns  About  a 
Bluff  Body  [Ref.  8] 


Trying  to  relate  the  flow  around  a  simple  bluff  body  to 
actual  airflow  over  ships  is  an  extremely  difficult  task. 
One  can  find  only  few,  and  then  poorly  done,  studies  on  the 
subject.  As  stated  before,  the  present  studies  at  NPS  should 
rectify  that  situation  soon.  Until  then,  it  is  necessary  to 
look  elsewhere.  One  such  area,  where  much  time  and  effort 
have  been  expended,  is  in  the  area  of  wind  flow  around 
buildings  and  other  obstacles  [Ref.  10]. 

Recent  investigations  in  the  flow  around  buildings  have 
advanced  the  understanding  of  physical  flow  processes 
occurring  in  the  near  and  far  wake  region.  But,  according  to 
Peterka  et  al.  [Ref.  10],  even  though  there  has  been  many 
studies  in  this  area,  there  is  still  a  high  level  of 
misunderstanding  on  how  the  winds  actually  flow  around 
buildings.  They  believe  that  these  misconceptions  are 
probably  caused  by  "conceptual  extensions  of  two-dimensional 
flow".  Figure  6  shows  the  separation  zone  for  a  two- 
dimensional  object  bounded  by  streamlines  so  that  the  cavities 
are  closed.  For  a  three-dimensional  object,  these  separation 
lines  are  no  longer  valid.  Studies  which  have  been  conducted 
in  boundary  layer  wind  tunnels  [Refs.  11,12]  show  that  three- 
dimensional  objects  show  fundamental  differences  in  flow 
patterns  relative  to  those  of  two-dimensional  objects.  Since 
this  study  is  indeed  based  on  two-dimensional  flow  relating  to 
a  three-dimensional  problem,  it  is  worthwhile  to  spend  a 
little  time  and  look  at  these  differences. 
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CORRECT  :  t-  0  BODIES 
INCORRECT :  3-D  BODIES 


Figure  6.  Separation  Cavities  for  2-d  Flow  [Ref.  10] 


Looking  at  Figures  7  and  8,  the  flow  approaching  the 
obstacle  has  separated  at  some  distance  upstream,  at  a  point 
that  is  dependent,  to  the  first  order,  on  building  height-to- 
width  ratio,  building  height-to-boundary-layer~height  ratio 
and  upstream  surface  roughness  (Ref.  10].  The  air  in  this 
separated  flow  strikes  the  building,  flows  downward  and  rolls 
up  into  a  vortex.  It  then  wraps  around  the  building  into  the 
horseshoe  shape  that  was  discussed  before.  This  horseshoe 
vortex  can  be  identified  in  the  flow  for  quite  some  distance 
downstream. 

The  wind  that  impinges  on  the  front  of  the  building  forms 
a  stagnation  region  somewhere  near  the  top  (about  2/3  the  way 
up)  depending  or.  building  height-to-width  ratio.  (Ref.  10) 
From  this  region,  flow  moves  out  toward  all  front  edges  of  the 
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Figure  7.  Centerline  Streamline  Patterns  for  Flow 
Reattaching  to  Top  [Ref.  10] 


Figure  8.  Centerline  Streamline  Patterns  for  Flow 
Not  Reattaching  to  Top  [Ref.  10] 


object.  Near  these  edges,  it  separates  and  may  or  may  not 
reattach  before  reaching  the  back  edge.  This  reattachment 
depends  on  many  factors  such  as  building-length-to  width 


ratio,  height-to-length  ratio  and  upstream  roughness  (which 
also  determines  the  turbulence  intensity  in  the  approaching 
wind) .  Figures  5  and  7  show  flow  patterns  for  attached  flow 
and  Figure  8  shows  flow  patterns  for  unreattached  flow. 

A  separation  cavity  covers  the  rear  face  of  the  object. 
The  cavity  length  is  defined  by  the  distance  from  the  building 
to  the  centerline  reattachment  point  downwind  of  the  body,  and 
is  normalized  by  the  building  height.  It  can  vary  from  two  to 
six  building  heights  [Ref.  8].  As  a  result  of  this  variation, 
the  reattachment  "point"  is  more  accurately  called  a  reattach¬ 
ment  zone.  This  subject  will  be  looked  at  more  closely  later. 

Obtaining  a  clear  picture  of  this  separation  cavity  is 
quite  difficult  due  to  the  high  level  of  turbulence  inside 
this  region.  Figure  9,  which  graphically  depicts  the  streak 
lines  of  helium  bubbles,  is  a  photograph  of  a  helicopter  deck 
from  a  flow-visualization  study  of  a  DD-963  class  destroyer 
model.  Bearing  in  mind  that  the  size  of  the  cavity  changes 
with  ship  yaw  angle  and,  to  some  degree,  with  pitching  and 
rolling  of  the  ship,  it  is  little  wonder  that  creating  an 
accurate  rotor  engage/disengage  envelope  has  been  both 
difficult  and  unsuccessful. 

C.  FLOW  OVER  A  BACKWARD  FACING  STEP 

The  next  logical  step  is  to  proceed  and  review  any  studies 
done  on  the  flow  over  a  backward  facing  step  (BFS) .  Not  only 
is  this  type  of  flow  probably  the  easiest  reattaching-flow 
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Figure  9.  Wake  of  Hangar;  DD-963  [Ref.  7:p.  73] 


scenario  to  observe  but,  it  also  comes  very  close  to 
resembling  the  flow  over  the  centerline  of  the  helicopter  deck 
on  the  back  of  the  AOR.  Luckily,  there  have  been  several 
different  type  studies  of  this  BFS  and  it  has  been  established 
as  a  "benchmark  problem"  for  checking  CFD  codes.  However,  the 
present  review  will  be  confined  to  those  which  were 
restricted  to  turbulent  two-dimensional  freestream  flow. 

Although  the  backward  facing  step  is  the  simplest  of  the 
reattaching  flows,  that  fact  is  in  no  way  indicative  of  the 
complexity  of  the  flowfield;  it  is  still  a  very  complex  flow, 
as  is  illustrated  by  Figure  10.  It  can  be  observed  that  the 
upstream  boundary  layer  separates  at  the  sharp  corner  forming 


a  shear  layer  [Ref.  13].  This  separated  shear  layer  curves 
sharply  downward  in  the  reattachment  zone.  Then  after 
striking  the  floor,  part  of  the  flow  is  deflected  upstream 
into  the  recirculating  flow  by  a  strong  adverse  pressure 
gradient.  It  would  be  incorrect  to  consider  this 
recirculation  area  as  a  dead  air  zone.  Backward  flow  has  been 
measured  at  over  20%  of  the  mean  freestream  velocity.  [Ref. 
13] 

Eaton  and  Johnston  [Ref.  13]  have  compared  the  results  of 
several  studies  [Refs.  14-16]  concerning  the  reattachment 
length,  which  together  with  the  level  of  turbulence  inside  the 
recirculation  area,  are  probably  the  most  important  parameters 
that  characterize  this  flowfield.  This  work  gives  insight 
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into  the  effect  of  varying  the  following  four  independent 
parameters:  initial  boundary-layer  state,  initial  boundary- 
layer  thickness,  freestream  turbulence  and  the  aspect  ratio. 

It  was  found  that  the  effect  of  changing  the  state 
(laminar/turbulent)  of  the  separation  boundary  layer  had  a 
significant  effect  on  the  reattachment  length  (see  Figure  11) . 
The  flow  apparently  becomes  independent  of  Reynolds  number, 
based  on  momentum  thickness,  when  the  boundary  layer  is  fully 
turbulent . 

Data,  in  Eaton  and  Johnston's  study,  show  that  the 
reattachment  length  has  a  weak  dependence  on  the  effect  of 
changing  the  state  (laminar/turbulent)  of  the  separation 
boundary  layer.  However,  four  other  data  sets  with  different 
values  of  the  boundary  layer  thickness,  but  similar  other 
parameters,  show  the  reattachment  length  having  a  much 
stronger  dependence  on  the  thickness.  The  data  suggest  that 
further  study  is  needed  to  resolve  this  issue. 

The  effect  of  freestream  turbulence  on  the  reattachment 
length  has  never  been  studied  systematically.  The  few  data 
sets  that  resulted  from  these  studies,  and  documented  in  Eaton 
and  Johnston's  study,  showed  that  fairly  high  levels  of 
turbulence  seemed  to  decrease  the  reattachment  length.  Again, 
further  investigation  is  required. 

The  effect  on  the  reattachment  length,  of  the  aspect  ratio 
of  the  flow  apparatus  (channel  height  to  step  height),  was 
also  documented.  These  studies  found  that  the  effect  was 
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Figure  11.  Reattachment-length  Measurements  Showing 
Dependence  on  the  State  of  the  Separating 
Boundary  Layer  JRef.  13] 


negligible  for  aspect  ratios  greater  than  ten.  For  aspect 
ratios  less  than  ten,  the  reattachment  length  increases  if  the 
boundary  layer  at  separation  is  laminar  and  decreases  if  it  is 


turbulent. 


The  other  important  parameter,  turbulence  in  the 
recirculation  area,  has  also  been  measured  for  most  of  the 
data  sets  in  Eaton  and  Johnston's  survey.  Though  there  seems 
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to  be  a  substantial  variation  in  the  peak  values  of  turbulence 
and  shear  stress,  the  turbulence-intensity  measurements  show 
a  consistent  pattern  when  the  maximum  intensity  is  plotted  as 
a  function  of  streamwise  distance.  In  almost  nil  cases,  the 
turbulence  intensity  reaches  its  greatest  value  approximately 
one  step  height  upstream  of  reattachment  and  then  decays 
rapidly  in  the  downstream  direction. 

D.  FLOW  OVER  A  TWO-DIMENSIONAL  OBSTACLE 

Since  a  deflector  will  be  used  to  direct  the  airflow  over 
the  backward  facing  step,  it  would  be  prudent  to  spend  a  brief 
time  looking  at  the  flow  over  a  two-dimensional  fence  immersed 
in  a  turbulent  layer  on  a  flat  surface.  Such  a  study  has  been 
conducted  recently  by  Atli  [Ref.  17]  who  analyzed  the  flow 
field  through  the  surface  oil  technique  of  flow  visualization. 
He  then  obtained  the  longitudinal  components  of  the  mean 
velocities  by  using  hot  wire  aneroometry  and  applied 
corrections  for  flow  reversal  and  turbulence. 

Figure  12  shows  the  structure  described  by  the  flow 
visualization  tests.  Primary  and  secondary  recirculation 
regions  exist,  both  upstream  and  downstream  of  the  fence,  with 
the  downstream  one  considerably  larger.  He  observed  that  the 
relative  height  of  the  obstacle,  in  terms  of  the  reference 
boundary  layer  and  consequently  the  Reynolds  number  based  on 
the  height  of  the  obstacle,  are  the  parameters  affecting  the 
shape  of  the  flowfield  and  the  structure  of  the  turbulence. 
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Figure  12.  Flow  Structure  Over  a  2-d  Fence  [Ref.  17] 


Specifically  the  data  indicated  that  the  relative  length  of 
the  primary  recirculation  region  slightly  increases  when  the 
relative  height  of  the  obstacle,  and  consequently  the  Reynolds 
number  based  on  obstacle  height,  increases.  He  did  not 
observe  this  effect  for  the  other  recirculation  areas;  in 
fact,  they  remained  the  same  size  for  the  changes  made. 

Atli  went  on  to  compare  his  work  with  that  of  Sinha  et  al. 
[Ref.  18]  who  worked  with  flow  over  a  backward  facing  step. 
He  discovered  that  for  2-d  flows  with  the  same  value  of 
Reynolds  number,  the  length  of  the  primary  recirculation 
region  downstream  of  a  fence  on  a  flat  surface  is  longer  than 
that  on  the  backward  facing  step.  The  reason  for  this  is  the 
flow  approaching  the  obstacle  diverges  from  the  horizontal 
before  reaching  the  obstacle  because  of  the  upstream 
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recirculation  region.  The  flow  approaching  the  backward 
facing  step  is  horizontal  and  does  not  diverge. 

Another  conclusion  which  Atli  reaches,  which  has  some 
relevance  to  the  present  study,  is  that  the  reverse  velocity 
profile  in  the  recirculation  region  increases  with  height  of 
the  step  and  therefore  with  the  height-based  Reynolds  number. 
This  increase  in  obstacle  height  and  Reynolds  number  also 
increases  the  maximum  turbulence  intensity  in  the  mixing 
region. 

E.  FLOW  MODIFICATION 

Though  no  studies  can  be  found  that  attempt  to  modify  the 
airflow  over  a  backward  facing  step  through  the  use  of 
deflectors,  there  have  been  several  studies  of  the  airflow 
over  and  through  fences  and  shelterbelts.  These  studies 
generally  deal  with  modifying  the  airflow,  so  as  to  protect 
crops  or  provide  comfort  for  humans,  through  the  use  of 
various  type  windbreaks.  Much  of  this  information  is 
applicable  to  our  current  study  and  will  be  used  in  the 
modelling  portion  of  the  problem. 

windbreaks  and  shelterbelts  have  played,  and  continue  to 
play,  an  important  part  in  protecting  man  and  his  environment. 
It  was  with  this  in  mind  that  the  World  Meteorological 
Organization,  at  its  second  session  in  1958,  first  set  up  a 
working  group  on  windbreaks.  They  worked  out  a  plan  for  long¬ 
term  experiments  to  assist  in  regional  planning  of  windbreaks 
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and  shelterbelts  for  research  purposes.  The  results  of  their 
work  were  published  in  1964  [Ref.  19]  and  much  of  what  they 
concluded  is  still  valid  today. 

One  of  the  general  conclusions  of  their  study  was  that 
beneath  the  peak  airflow  over  a  windbreak  is  the  zone  of 
greatest  wind  reduction  (see  Figure  13).  At  that  time  it  was 
called  a  "dead  calm  area"  by  Kreutz.  [Ref.  19 :p.  72]  This 
term  is  really  a  misnomer,  as  has  since  been  shown.  [Ref.  17] 
Figure  14  illustrates  the  reduction  in  horizontal  wind  behind 
a  windbreak  as  a  function  of  windbreak  permeability.  This 
figure  also  shows  that  the  lower  the  porosity  of  the  obstacle, 
the  nearer  to  the  obstacle  the  "calm"  area  is  located. 

The  crucial  parameter  for  wind  reduction  behind  shelter- 
belts  was,  therefore,  determined  to  be  the  shelterbelt *s 
density,  or  porosity.  The  less  porous,  the  more  wind 
reduction,  but  for  only  a  very  small  area  immediately  behind 
the  shelterbelt.  Then,  as  the  porosity  increases,  the 
horizontal  wind  velocity  increases  slightly,  but  the  area  of 
protection,  measured  downwind,  also  increases.  The  overall 
best  protection,  which  extended  six  to  seven  shelterbelt 
heights  downstream,  was  thought  to  be  about  50%  porosity. 

In  1981  Perera  [Ref.  20]  showed  that  the  normalized  mean 
wind  velocity  through  a  porous  fence  was  independent  of  the 
form  of  the  fence  construction.  Small  holes,  large  holes,  and 
even  horizontal  slat  fences  seemed  to  have  no  significant 
effect  on  altering  the  results  for  a  given  porosity,  several 
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Figure  13.  Mean  Streamlines  at  a  Medium  Dense  and  a 
Dense  Windbreak  [Ref.  19 :p.  120] 


other  studies  have  also  confirmed  Perera's  results.  Another 
of  his  conclusions  was  that,  as  the  porosity  increased,  the 
recirculation  bubble  decreased  in  size  and  moved  downstream. 
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His  results  showed  that  the  recirculation  bubble  existed  only 
for  fences  with  porosities  less  than  30%. 


Figure  14.  The  wind  Speed  Reduction  of  Different  Wind¬ 
breaks  (H  =  Height  of  Windbreak)  [Ref.  19: 

p.  120] 


To  summarize,  wake  velocities  increase  with  increasing 
porosity  and,  in  contrast,  the  turbulence  intensity  decreases 
with  increasing  porosity. 

Perera  uses  a  parameter  suggested  by  Gandemer  [Ref.  21], 
which  takes  into  account  both  the  mean  and  turbulence 
properties  of  the  wind  to  determine  the  best  overall  porosity 
for  a  fence.  He  (Perera)  suggests  a  fence  with  a  porosity  of 
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about  10%.  Interestingly,  Gandemer  cc  icluded  that  a  fence  of 
20%  porosity  provided  the  best  protection,  but  that  a  fence 
of  50%  porosity  was  almost  as  effective  and  less  expensive  to 
build.  It  is  hoped  that  all  of  these  porosity  parameters  can 
be  compared  during  the  computer-modeling  process. 
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III.  CFD  PROGRAM  SETUP 


A.  BACKGROUND 

The  computational  Fluid  Dynamics  code  used  for  the 
modelling  process  was  the  PHOENICS  Flow-simulation  code.  This 
code  was  first  created  in  1979-80  and  released  in  1981.  [Ref. 
22 :p.  l.i]  it  is  a  computer  code  which  simulates  fluid  flow, 
heat  transfer,  chemical  reactions  and  related  phenomena.  The 
latest  version  of  PHOENICS,  1.5,  dated  September  1989  was 
utilized  as  well  as  PHOENICS  "Easy  Flow,"  a  problem  definition 
and  management  program.  Easy  Flow  uses  a  menu-driven, 
graphical  interface  for  problem  setup,  which  is  interfaced 
with  the  PHOENICS  code.  This  greatly  simplified  setting  up 
the  problems  by  eliminating  the  tedious  and  long  task  of 
writing  new  code  each  time  a  minor  change  was  made  to  the 
problem. 

PHOENICS  consists  of  two  essential  computer  codes  and  four 
auxiliary  ones.  The  essential  codes  consist  of  a  pre¬ 
processor  called  SATELLITE  and  a  processor  called  EARTH.  The 
auxiliary  codes  are  post-processors  called  PHOTON,  AUTOPLOT 
and  PINTO,  and  a  separate  self-instruction  program  called 
GUIDE.  [Ref.  22:p.  2.3] 

SATELLITE  is  an  interpreter;  from  instructions  provided  by 
the  user  it  creates  a  data  file,  Q1  (see  Appendix) ,  containing 
instructions  which  EARTH  can  obey. 
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EARTH  contains  the  main  flow-simulating  software.  It 
incorporates  coding  sequences  which  represent  the  relevant 
laws  of  physics  applied  to  elements  of  material  distributed  in 
space  and  time  [Ref.  22:p.  2.3].  EARTH  reads  the  Q1  file 
provided  by  SATELLITE  and  executes  the  corresponding 
calculations.  It  then  produces  an  output  file  called  RESULT, 
which  the  user  can  read  or  have  sent  to  another  file  called 
PHIDA. 

The  PHOTON  code  is  an  interactive  program  which  picks  up 
the  PHIDA  file  written  by  EARTH  and  then,  in  response  to 
instructions  entered  by  the  user  through  the  VDU  keyboard, 
represents  the  computed  grid  and  flow  pattern  graphically  on 
the  screen. 

B.  HOW  "PHOENICS"  WORKS 

PHOENICS  describes  phenomena  involving  the  flow  of  heat 
or  material  in  terms  of  distributions  in  space  and  time  of 
temperatures,  velocities,  pressures,  and  other  physical 
quantities.  These  are  the  dependent  variables.  These 
distributions  involve  ascribing  numerical  values  to  the 
variables  at  each  of  an  orderly  array  of  locations,  called 
"nodes"  or  "grid  points."  [Ref.  22 :p.  2.4] 

The  printed  or  plotted  values  of  temperatures,  velocities, 
etc.,  produced  by  PHOENICS  represent  solutions  of  a  set  of 
algebraic  equations  representing  the  relevant  laws  of  physics 
for  the  chosen  set  of  cells  (grid  points)  (see  Figure  15)  . 
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These  laws  are  those  governing  the  conservation  of  mass, 
momentum  and  energy,  coupled  with  the  transport  laws 
describing  heat  conduction,  diffusion  of  matter,  and  viscous 
action  [Ref.  22:p.  2.4]. 


Figure  15.  Cell  Depiction  [Ref.  22:p.  1.11] 


In  the  present  study,  no  conduction  nor  diffusion  nor 
reactions  are  considered  and  the  relevant  equations  are  the 
incompressible,  constant  viscosity  time  averaged  Navier-Stokes 
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equations;  discretized  via  a  control  volume  formulation  into 
the  generic  form; 


Ve  +  Vw  +  aN*N  +  Vs  VH  +  aL»L  *  Vt  +  S  , 

P  *  aE  +  aW  +  »N  +  aS  +  aH  +  aL  +  aT  +  ®P 

where  tp  represents  the  dependent  variable  currently  in 
question.  The  subscripts  P,  E,  W,  N,  H,  and  L  denote  the 
location  at  which  this  variable  is  computed,  in  accordance 
with  the  north-south-east-west-high-low  convention  used  in 
PHOENICS  [Ref.  22 :p.  2.2]. 

For  each  dependent  variable  there  are  as  many  algebraic 
equations  as  there  are  cells  in  the  discretized  grid  domain. 
For  example,  in  a  typical  two-dimensional  problem  with  100 
cells  in  both  the  X  and  Y  directions,  and  five  dependent 
variables  (from  conservation  equations),  there  would  be  50,000 
algebraic  equations  to  be  solved.  In  addition  to  being 
numerous,  these  sets  of  equations  are  strongly  coupled  through 
convective  terms.  PHOENICS  solves  the  discretized  equations 
through  an  iterative  method  reducing  the  imbalance  between  the 
right  and  left  sides  of  every  equation  to  a  magnitude  which  is 
small  enough  to  be  considered  negligible  (i.e.,  a  converged 
solution).  [Ref.  22:p.  2.13] 

The  iterative  process  itself  is  very  complex,  involving  a 
multi-stage  sequence  of  adjustments  of  values,  repeated  many 
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times.  For  the  problem  in  this  study  a  "slabwise"  solution 
was  directed.  "Slabs"  are  arrays  of  cells  having  the  same 
value  of  the  low-to-high  coordinate  z.  A  "sweep"  is  a  set  of 
slab-wise  operations  conducted  in  sequence  from  the  lowest 
slab  to  the  highest.  Because  the  equations  for  values  at  one 
slab  ordinarily  make  reference  to  values  at  the  next  higher 
slab,  later  adjustments  made  at  the  higher  slab  will 
invalidate,  to  some  extent,  the  adjustments  which  have  just 
been  made  at  the  lower  one  [Ref.  22  :p.  2.11].  It  becomes 
obvious  then,  that  to  get  a  reasonable  and  accurate  solution, 
many  sweeps  must  be  made  until,  ideally,  all  equations  are  in 
a  "perfect  balance"  (i.e.,  a  converged  solution)  that  further 
adjustments  are  unnecessary. 

C.  PROBLEM  CONFIGURATION 

The  goal  of  this  paper  is  to  recommend  possible 
configurations  to  be  used  in  the  next  phase  of  the  program, 
which  is  an  experimental  verification  of  the  present  results, 
using  flow  visualization  and  measurement  techniques,  in  the 
NPS  low  speed  wind  tunnel.  To  that  end,  the  problem 
dimensions  used  in  PHOENICS  were  tailored  as  closely  as 
possible  to  the  equipment  and  physical  dimensions  available  in 
the  wind  tunnel. 

The  backward  facing  step  has  a  height  of  125mm  and  the 
deflectors  ranged  in  length  from  25  to  50mm.  All  of  these 
dimensions  were  incorporated  into  the  PHOENICS  program.  The 
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freestream  velocity  and  turbulence  profiles  -ere  measured  in 
the  wind  tunnel .  These  velocity  and  turbule  ice  profiles  are 
given  in  Table  1  and  represent  approximately  a  20  m/sec  wind 
over  a  fully  developed  sea. 


TABLE  1 

FLOW  SIMULATION  INPUTS 


height  (m) 

velocity  (m/sec) 

turbulent  intensity  (%) 

0.050 

1.91 

11.85 

0.075 

1.92 

12.32 

0.100 

2.05 

11.33 

0.200 

2.36 

4.79 

0.300 

2.40 

3.92 

0.400 

2.47 

3.66 

0.475 

2.56 

3.24 

0.625 

2.62 

3.18 

0.750 

2.72 

2.19 

Freestream*  kinetic  energy  -  k.  •  2.58  m*/s* 

Freestream*  velocity  -  u„  ■  2.77  m/sec 
Freestream*  dissipation  rate  -  c„  -  .172  mV®« cJ 
Boundary  layer  *  £  -  .762  m 

*  value  at  top  of  boundary  layer 

The  turbulence  model  selected  is  the  k ~c  two  equation  one. 

This  model  has  been  used  extensively  and  the  results  of 

Murakami  and  Mcchida  [Ref.  24]  show  that  it  yields  a 

reasonable  representation  of  the  velocity  field  around  a  cubic 
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model.  A  problem  with  the  k-£  model  was  revealed  by  an 
experiment  conducted  by  Yeung  and  Scott.  [Ref.  25:p.  27]  In 
that  study,  the  reattachment  lengths  obtained  were  under¬ 
predicted  when  compared  to  actual  case  studies.  They  believe 
this  to  be  an  inherent  weakness  of  the  two  equation  model. 

The  turbulence  viscosity  is  determined  from  the  values  of 
the  turbulent  kinetic  energy  k  and  its  dissipation  rate  c 
according  to: 


w 


t 


p  C, 


e 


where  p  is  the  fluid  density  and  is  a  constant. 

The  values  of  k  and  r  are  determined  from  their  own 
transport  equations: 
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The  values  of  CM,  C, ,  C,,  and  the  other  parameters  are  given  in 
the  nomenclature  list. 
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To  establish  a  baseline  on  how  the  results  of  PHOENICS 
compare  to  other  2-d  flow  studies,  two  basic  configurations 
were  used.  The  first  was  a  backward  facing  step  without  a 
deflector,  for  which  the  inputs  are  detailed  in  Table  1. 
Next,  solid  and  porous • fences  were  considered,  the  program 
run,  and  the  results  recorded.  The  porosity,  which  is  the 
percentage  of  the  deflector  area  open  to  airflow,  ranged 
between  15%  and  40%. 

Deflectors  were  then  introduced  into  the  study  in  an 
attempt  to  tailor  the  airflow  in  such  a  way  as  to  reduce  the 
velocity  gradient  and  the  overall  turbulence  in  the  region 
downstream  of  the  step.  The  first  such  deflector  was  50mm  in 
length  and  placed  25mm  directly  above  the  step's  edge  (see 
Figure  16)  .  This  configuration,  where  configuration  refers  to 
size,  location  and  orientation  of  a  deflector,  was  run  at 
different  porosities  before  a  new  configuration  was 
established  and  a  set  of  new  runs  attempted.  All  of  the  runs, 
including  deflector  size,  position,  and  porosity  are 
documented  in  Table  2 . 

Next,  in  order  to  create  a  slanted  deflector.  Body  Fitted 
Coordinates  (BFC)  were  utilized.  An  example  of  such  a 
configuration  is  shown  in  Figure  17.  The  deflector  was  set  at 
angles  of  30  degrees  and  4  5  degrees  (measured  from  the 
vertical)  at  different  locations  and  porosities  for  several 
runs  (see  Table  2)  . 
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TABLE  2 


DEFLECTOR  CONFIGURATIONS 


Distance  Offset  Dis- 
Above-Below  tance  Down- 

Deflector/  Porosity  Step  Edge*  stream  of  Step 


Figure  Fence-length 

(%) 

(mm) 

Face*  (mm) 

18a 

18b 

- 

- 

- 

- 

19a 

19b 

125 

0.0 

- 

- 

20a 

20b 

125 

10.0 

- 

- 

2la 

21b 

125 

20.0 

- 

- 

22a 

22b 

50 

15.0 

25.0 

0.0 

23a 

23b 

50 

0.0 

12.5 

12.5 

24a 

24b 

50 

15.0 

12.5 

12.5 

25a 

25b 

50 

0.0 

0.0 

12.5 

26a 

26b 

50 

5.0 

12.5 

25.0 

27a 

27b 

25 

0.0 

0.0 

50.0 

28a 

28b 

50 

0.0 

0.0 

50.0 

29a 

29b 

50 

5.0 

0.0 

50.0 

30a 

30b 

50 

15.0 

-25.0 

50.0 

31a 

31b 

50 

0.0 

0.0 

30.0 

32a 

32b 

50 

10.0 

0.0 

30.0 

33a 

33b 

50 

0.0 

0.0 

50.0 

34a 

34b 

50 

5.0 

0.0 

50.0 

35a 

35b** 

50 

0.0 

0.0 

50.0 

36a 

36b*** 

50 

0.0 

0 . 0 

50.0 

* 

Distance 

measured 

from  lower 

edge  of 

deflector 

**  30/60  degree  cranked  deflector 

***  45/60  degree  cranked  deflector 
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IV.  RESULTS  AND  DISCUSSION 


The  configurations  used  in  the  "Phoenics"  CFD  program  are 
documented  in  Table  2.  Briefly,  they  were:  a  backward  facing 
step;  a  fence  set  perpendicular  to  the  flow  and  various 
deflector  arrangements  in  conjunction  with  the  backward  facing 
step.  All  runs  were  allowed  to  continue  the  iteration  process 
until  such  time  that  the  residuals  were  small  enough,  or  had 
decreased  in  magnitude  enough,  to  achieve  what  was  considered 
balanced  equations  (converged  solution).  [Ref.  25  p.  5.26] 
Residuals  refer  to  the  measure  of  total  imbalance  of  the 
discretized  equations  being  computed.  The  number  of  sweeps  to 
accomplish  this  varied  between  600  and  900. 

Initially  the  code  as  installed  on  the  COMPAQ  DESKPRO 
386/25,  would  complete  only  a  small  number  of  sweeps  (between 
one  and  100)  before  an  error  statement  appeared  on  the  VDU  and 
all  iterations  ceased.  After  several  weeks  of  troubleshooting 
and  fruitless  calls  for  assistance  to  CHAM  (the  developers  of 
"Phoenics") ,  it  was  decided  to  remove  the  WEITEK  math  co¬ 
processor  which  was  installed  in  the  computer  along  with  the 
Intel  80387  math  co-processor.  CHAM  made  assurances  that  the 
software  package  sent  was  specifically  designed  to  be  used 
with  both  co-processors  and  that  removing  the  WEITEK  co¬ 
processor  would  have  no  effect.  However,  when  the  removal  was 
complete,  the  "Phoenics''  code  functioned  normally.  Other 
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problems  encountered  with  the  "Phoenics"  code  will  be 
discussed  later. 

A.  BACKWARD  FACING  STEP 

The  first  runs  were  made  simulating  airflow  over  a  125mm 
high  backward  facing  step  (BFS)  as  shown  in  Figure  18a  and 
18b.  The  velocity  vectors  and  contours  of  turbulent  kinetic 
energy  (TKE)  values  were  plotted  for  all  configurations. 

The  flow  over  the  BFS  closely  resembles  that  which  has 
been  documented  in  the  literature  previously  referenced,  with 
one  exception,  the  reattachment  zone.  It  is  located  only 
about  3.5  step  heights  downstream  rather  than  the  four  to  six 
step  heights  achieved  in  wind  tunnel  experiments.  This  is 
more  than  likely  a  result  from  using  the  k-e  turbulence  model 
which,  as  was  stated  previously,  tends  to  underestimate  the 
reattachment  zone.  However,  on  a  positive  note,  the  area  of 
maximum  TKE  is  located  approximately  one  step  height  upstream 
from  the  reattachment  zone,  as  was  the  case  in  the  literature 
cited . 

These  velocity  vectors  and  TKE  values  of  Figure  18  will 
now  be  used  as  a  '•benchmark"  for  comparison  against  all  other 
configurations  used.  A  region  of  special  interest  to  this 
study  is  the  region  of  space  that  would  be  occupied  by  the 
main  rotor  blades  of  the  helicopters  during  start-up  and  shut¬ 
down.  Hereafter  this  region  will  be  referred  to  as  "the 
region  of  special  interest."  It  ranges  from  one-half,  to  one 


39 


full  step  height  downstream  of  the  step  and  one-fourth,  to 
three-fourths  of  a  step  height  off  the  floor.  In  the 
discussion  that  follows,  unless  otherwise  indicated, 
comparisons  of  the  properties  of  a  given  flow  and  the 
reference  BFS  flow  refer  to  this  particular  area. 

B.  TWO-DIMENSIONAL  FENCE 

The  flow  depicted  downstream  of  the  fence  has  a  well- 
defined  recirculation  area  with  reverse  flow  up  to  30%  greater 
than  that  found  downstream  of  the  BFS,  and  TKE  intensities  up 
to  400%  greater,  as  shown  in  Figure  19a  and  19b. 

Runs  with  a  10%  porous  fence  resulted  in  the  primary 
recirculation  area  shifting  approximately  1/2  step  heights 
downstream  and  a  secondary  recirculation  zone  appearing 
immediately  downstream  of  the  fence,  as  shown  in  Figure  20a 
and  20b.  Evidence  of  this  secondary  zone  downstream  of  a 
porous  fence  has  not  been  found  in  the  references  cited.  The 
TKE  intensities  increased  greatly  in  the  vicinity  of  the  fence 
and  then  quickly  dissipated  as  the  flow  moved  downstream. 

As  shown  by  Figure  21a  and  21b  the  20%  porous  fence 
completely  eliminated  all  traces  of  the  recirculation  area  but 
the  fence  had  even  higher  values  of  TKE  around  it.  Most 
literature  cited  had  the  recirculation  area  existing  until  at 
least  30%  porosity  was  achieved.  This  discrepancy  will  have 
to  be  further  examined  when  the  study  moves  to  the 
experimental  stage. 
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Figure  20. 


2-d  Fence;  10%  Porosity 


C.  VERTICAL  DEFLECTOR  ABOVE  STEP 


The  first  configuration  using  a  vertical  deflector  was  one 
in  which  a  50mm  deflector  was  placed  25mm  directly  above  the 
step.  At  0%  porosity  there  was  a  great  deal  of  recirculation 
behind  the  deflector,  while  at  15%  porosity  this  recirculation 
was  eliminated,  as  illustrated  in  Figure  22a  and  22b.  In  both 
cases,  there  was  little  c"  no  effect  on  the  recirculation 
behind  the  step.  Decreasing  the  vertical  gap  (measured  from 
the  top  edge  of  the  step  to  the  bottom  edge  of  the  deflector) 
to  12.5mm  had  the  same  overall  results. 

D.  VERTICAL  DEFLECTOR  OFFSET  12 . 5MM 

With  the  vertical  gap  still  at  12.5mm  the  50mm  deflector 
was  offset  downstream  of  the  step  by  12.5mm.  Figure  23a  and 
23b  indicate  that  at  0%  porosity,  a  recirculation  area  exists 
that  is  much  larger  than  that  found  behind  the  BFS ,  though  it 
was  shifted  downstream  approximately  one  additional  step 
height.  The  velocities  in  the  "region  of  interest"  were,  on 
the  average,  of  the  same  magnitude  as  those  of  the  BFS,  though 
the  flow  vectors  were  now  all  moving  in  the  same  upward 
direction.  The  TKE  was  very  high  in  the  vicinity  of  the 
deflector  but,  near  the  floor,  was  lower  than  for  the  BFS. 

Increasing  the  porosity  to  5%,  10%  and  then  15% 
increasingly  reduced  the  recirculation  area.  At  15%  porosity, 
the  primary  recirculation  area  resembles  that  of  the  BFS? 
however,  the  velocity  vectors  were  all  reduced  10%-20%  as 
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Figure  22.  50mm  Vertical  Deflector;  15%  Poroeity 
25mm  Vertical  Gap 


shown  in  Figure  24a  and  24b.  The  TKE  was  also  reduced  in  "the 
region  of  interest." 

Lowering  the  deflector  so  that  there  was  no  vertical  gap 
resulted  in  a  very  large  recirculation  area  with  velocity 
vectors  20%-30%  larger  than  those  found  downstream  of  the  BPS 
(see  Figure  25) .  Interestingly,  the  TKE  intensities  remained 
about  the  same  as  the  BFS . 

It  became  obvious  that,  to  reduce  the  primary  recircula¬ 
tion  zone,  the  deflector  would  have  to  be  moved  further 
downstream. 

E.  VERTICAL  DEFLECTOR  OFFSET  25MM 

Moving  the  50mm  deflector  to  a  position  25mm  downstream 
with  a  12.5mm  vertical  separation  and  0%  porosity  did  not  give 
much  better  results.  However,  increasing  the  porosity  to  5%, 
eliminated  most  of  the  recirculation  downstream  of  the 
deflector  and  reduced  the  velocity  vectors  downstream  of  the 
step  on  the  average  of  about  70%,  as  illustrated  in  Figure  26a 
and  26b.  This  drastically  reduced  the  velocities  in  the 
recirculation  zone.  While  the  TKE  is  quite  high  in  the 
vicinity  of  the  deflector,  it  dissipates  rapidly  as  the 
distance  increases,  and  becomes  much  lower  than  that  of  the 
BFS  in  "the  region  of  interest." 
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F.  VERTICAL  DEFLECTOR  OFFSET  50MM 

The  deflector  was  moved  further  downstream  to  a  position 
50mm  from  the  step  and  without  any  vertical  gap.  Both  25mm 
and  50mm  deflectors  were  used  with  different  results. 

The  25mm  deflector  ^t  0%  porosity  channeled  the  airflow 
in  such  a  way  as  to  eliminate  the  primary  recirculation  zone, 
as  seen  in  Figure  27a  and  27b,  but  a  smaller  and  weaker 
recirculation  zone,  with  little  TKE  was  formed  at  tne  base  of 
the  step.  In  "the  region  of  interest"  the  flow  is 
predominantly  downward  with  a  velocity  magnitude  such  that  is 
found  at  the  downstream  end  of  the  BFS  recirculation  zone. 
The  TKE  is  approximately  the  same  as  that  found  in  the  BFS, 
but  was  much  lower  near  the  floor.  Increasing  the  porosity, 
even  to  5%,  increased  the  airflow  enough  to  reform  the  large 
primary  recirculation  zone. 

Increasing  the  size  of  the  deflector  to  50mm  and  setting 
the  porosity  to  0%  eliminated  the  secondary  recirculation  zone 
at  the  base  of  the  step,  but  another  small  zone  appeared 
slightly  below  and  downstream  of  the  deflector;  this  is 
illustrated  by  Figure  28a  and  28b.  The  flow  is  now  directed 
downward  with  localized  velocities  under  the  deflector  and 
with  the  TKE  up  to  four  times  that  found  downstream  of  the 
BFS. 

By  increasing  the  porosity  to  5%  the  secondary  recircula¬ 
tion  zone  is  eliminated,  but  there  still  remains  a  low 
velocity  "S"  type  flow  pattern,  as  seen  in  Figure  29a  and  29b. 
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The  TKE  is  slightly  lower  than  the  previous  example  but  still 
far  greater  than  what  is  foi’d  downstream  of  the  BFS . 
Increasing  the  porosity  any  further  allows  the  primary 
recirculation  zone  to  reappear  with  both  velocity  vectors  and 
TKE  substantially  higher  than  the  BFS. 

An  additional  configuration  was  considered,  in  which  a 
50mm  deflector  was  again  offset  50mm,  but  now  placed  such  that 
half  of  the  deflector  was  above  the  step  and  1/2  below.  The 
only  run  that  had  an  interesting  effect  was  the  15%  porous 
deflector.  The  primary  recirculation  zone,  as  shown  in  Figure 
30a  and  30b  was  almost  eliminated,  though  a  small,  weak, 
remnant  is  still  evident  further  downstream.  An  "S"-type  flow 
pattern  is  formed  with  velocity  and  TKE  magnitudes 
approximately  the  same  as  the  BFS.  However,  the  velocity 
vectors  do  not  start  their  upward  flow  until  much  further 
downstream. 

G.  SLANTED  OFFSET  DEFLECTORS 

The  first  such  slanted  deflector  was  offset  30mm  and 
angled  back  30  degrees  (see  Figure  31)  .  Remnants  of  the 
original  primary  recirculation  zone  are  just  visible  at  the 
right  edge  of  the  figure.  Another  "S"-type  flow  can  be 
observed  with  velocities  of  about  the  same  magnitude  as  those 
of  the  BFS;  in  essence,  there  are  two  recirculation  zones  with 
opposite  sense,  one  being  stronger  than  the  other.  The  TKE  in 
"the  region  cf  interest"  is  much  lower  than  that  of  the  BFS. 
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Figure  30.  50mm  Vertical  Deflector;  15%  Porosity 
Bottom  Edge  -25mm  Below  Step  Top, 

50mm  Offset  Downstream 


It  is  interesting  to  note  that  "Phoenics"  does  not  predict 
any  measurable  flow  immediately  along  the  back  of  the  step 
while  using  the  angled  deflectors.  This  was  not  the  case  when 
using  the  straight  vertical  deflectors. 

Increasing  the  porosity  to  10%  and  above,  increased  the 
size  of  the  displaced  primary  recirculation  zone  so  that  there 
are  now  two  significant  zones  of  recirculation  to  contend 
with.  These  are  illustrated  in  Figure  32a  and  32b.  The  TKE, 
apart  from  the  vicinity  of  the  deflector,  is  much  lower  than 
for  the  BFS . 

The  offset  distance  was  increased  to  50mm,  and  at  0% 
porosity  tl  primary  recirculation  zone  is  again  shifted 
downstrea  j  shown  in  Figure  33a  and  33b.  It  appears  to  be 
extremely  weak  and  of  little  consequence.  The  "S"  pattern  is 
again  evident  but  observing  "the  area  of  interest"  the  flow  is 
fairly  uniform  with  the  TKE  approximately  the  same  magnitude 
as  the  BFS. 

On  increasing  the  porosity  to  5%,  not  only  has  the  primary 
recirculation  zone  become  more  pronounced,  but  the  "S"  flow 
velocity  vectors  have  increased  on  the  order  of  25%-40%  with 
a  corresponding  increase  in  the  TKE  (see  Figure  34)  .  A 
similar  trend  occurs  at  higher  porosities. 

The  final  two  configurations  were  cranked-type  deflectors. 
In  the  first  instance  the  lower  part  of  he  deflector  is  angled 
at  30  degrees  to  the  vertical  and  the  upper  at  60  degrees  to 
the  vertical,  as  shown  in  Figure  35a  and  35b.  The  second  case 
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differed  from  the  first  in  that  the  lower  deflector  was  angled 
at  45  degrees  to  the  vertical,  as  shown  in  Figure  36a  and  36b. 

Both  configurations,  however,  gave  quite  similar  results. 
Both  had  the  primary  recirculation  zones  shifted  downstream 
with  a  strong  "S"  type  flow  located  in  the  center  of  "the  area 
of  interest."  Again,  the  TKE  remained  at  approximately  the 
same  magnitude  as  found  in  the  BFS . 

Different  porosities  were  tried  but  the  iterations  were 
interrupted  each  time  by  a  numerical  overflow.  It  was  obvious 
that  the  automatic  relaxation  features  built  into  "Phoenics" 
were  not  capable  of  dealing  with  this  level  of  complexity. 
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V.  CONCLUSIONS-  AND  RECOMMENDATIONS 


The  results  of  this  study  show  that  a  50mm,  5%  porous 

vertical  deflector,  whose  lower  edge  is  25mm  downstream  of, 
and  12.5mm  above  the  upper  corner  of  the  step,  provided  the 
best  overall  reduction  of  those  deflectors  tested,  in  velocity 
and  turbulent  kinetic  energy  for  a  125mm  high  backward  facing 
step.  The  velocity  has  been  reduced  on  an  average  of  50%  to 
70%  and  the  turbulent  kinetic  energy  has  been  reduced  upwards 
to  300%. 

The  porosity  of  the  deflectors  appears  to  be  an 
independent  parameter  and  therefore  no  one  recommendation  can 
be  made.  Each  deflector  position  appears  to  require  a  unique 
porosity  to  achieve  the  best  possible  airflow. 

It  is  also  clear  that  incorrect  deflector  placement  or 
orientation  can  do  more  harm  than  good  by  actually  increasing 
both  the  velocities  and  the  THE. 

Cranked  deflectors,  made  of  two  attached  straight 
segments,  appear  to  be  no  improvement  over  single  slanted 
segments . 

Some  recommendations  for  further  studies  are: 

(1)  Attempt  using  as  fine  a  grid  as  possible  when  setting 
up  problems.  This  will  take  more  computer  time  but 
give  a  more  detailed  and  accurate  representation  of  the 
flowf ield. 
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(2)  Rerun  the  program  for  lower  freestream  velocities  and 
compare  the  results  with  this  present  study's  results 
to  determine  any  relationships. 

(3)  Use  curvilinear  Body  Fitted  Coordinates  in  "Phoenics" 
to  set  up  curved  type  deflectors  to  learn  more  on  the 
effect  of  deflector  shape  on  airflow  behavior. 

(4)  Extend  the  study  to  3 -dimensional  flow,  i.e.,  study  the 
influence  of  the  airflow  coming  over  the  helicopter 
deck  from  around  the  sides  of  the  hangar  and  especially 
from  over  the  sides  of  the  ship. 

Complex  flow  problems  that  result  in  numerical  overflow 
errors  will  require  adjusting  the  "relaxation  factors”  in  the 
"Phoenics'1  code  so  that  convergence  can  occur. 
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APPENDIX 


SAMPLE  fiLEILE 


TALK=T ; RUN (1,1); VDU= 1 

♦A********************************************************** 

* 

*  GROUP  l.  Run  title  and  other  preliminaries. 

* 

TEXT (  Slanted  Deflector  ) 

* 

★A********************************************************** 

* 

*  GROUP  2.  Transience;  time-step  specification. 

* 

* 

* 

*  GROUP  3.  X-direction  grid  specification. 

* 

*  Cartesian  Grid  Selected 
CARTES=T 

*  Extent  of  the  Domain  in  the  X-Direction:  1.276E+00 

*  Number  of  Cells  in  the  X-Direction 
NX=29 

*  Equal  Grid  Spacing  in  Subregion  1 

*  SUBGRD(X, 1,8,  3 . 520E-01 ,  1.000E+00) 
INTEGER(NXF01,NXL01) ;  NXF01=1;  NXL01=8 
XFRAC ( 1)  =  4 . 4  00E-02 ; XFRAC ( 2 )  =  8.800E-02 
XFRAC(3)=  1. 320E-01;XFRAC(4)=  1.760E-01 
XFRAC (5 ) =  2 . 200E-01 ; XFRAC ( 6) =  2.640E-01 
XrRAC(7)=  3.080E-01; XFRAC ( 8 )  =  3.520E-01 

*  Equal  Grid  Spacing  in  Subregion  2 

*  SUBGRD  ( X ,  9 , 10 ,  8.800E-02,  1.000E+00) 

INTEGER ( NXF02, NXD02  )  ;  NXF02=9 ;  NXL02=10 
XFRAC  (9)  *=  3.960E-01;  XFRAC  ( 10)  =  4.400E-01 

*  Equal  Grid  Spacing  in  Subregion  3 

*  STJBGRD(X,  11, 11,  4 . 400E-02  ,  1.000E+00) 

INTEGER (NXF03 ,NXL03) ;  NXF03=11;  NXL03=11 
XFRAC (11)=  4 . 840E-01 

*  Equal  Grid  Spacing  in  Subregion  4 

*  SUBGRD(X, 12,29,  7.920E-01,  1.000E+00) 

INTEGER (NXF04, NXL04 ) ?  NXF04=12?  NXL04=29 
XFRAC (12)*=  5 . 280E-01 ; XFRAC ( 13 )  =  5.720E-01 
XFRAC  ( 14  )  ==  6 . 160E-01 ;  XFRAC  (15)  =  6.600E-01 
XFRAC (16)=  7.040E-01 ; XFRAC ( 17 ) =  7.480E-01 
XFRAC (18)=  7.920E-01 ; XFRAC ( 19 ) =  8.360E-01 
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XFRAC (20)  =  8.800E-01; XFRAC (21)=  9.240E-01 
XFRAC  (22)  =  9 . 68 OE-O 1 ; XFRAC  ( 23 )  =  1.012E+00 
XFRAC ( 24 )  =  1. 056E+00;XFRAC(25)=  1.100E+00 
XFRAC ( 26)  =  1. 144E+00;XFRAC(27)  =  1.188E+00 
XFRAC (28 )  =  1.232E+00;XFRAC(29)=  1.276E+00 

* 

************************************************************ 

* 

*  GROUP  4.  Y-direction  grid  specification. 

* 

*  Extent  of  the  Domain  in  the  Y-Direction:  1.276E+00 

*  Number  of  Cells  in  the  Y-Direction 
NY=29 

*  Equal  Grid  Spacing  in  Subregion  1 

*  SUBGRD(Y, 1,3,  1.320E-01,  1.000E+00) 

INTEGER (NYF01,NYL01) ;  NYF01=1;  NYL01=3 
YFRAC ( 1 ) =  4 . 4  OOE-02 ; YFRAC ( 2 )  =  8.800E-02 
YFRAC ( 3 )  =  1 . 320E-01 

*  Equal  Grid  Spacing  in  Subregion  2 

*  SUBGRD ( Y , 4 , 4 ,  4.400E-02,  1.000E+00) 

INTEGER  (NYF02,  NYL02  )  ;  NYF02=4 ;  NYL02=4 
YFRAC (4 ) =  1.760E-01 

*  Equal  Grid  Spacing  in  Subregion  3 

*  SUBGRD ( Y , 5 , 5 ,  4.400E-02,  1.000E+00) 

INTEGER (NYF03, NYL03 ) ?  NYF03=5;  NYL03=5 
YFRAC ( 5 ) -  2 . 200E-01 

*  Equal  Grid  spacing  in  Subregion  4 

*  SUBGRD (Y ,6,6,  4.400E-02,  1.000E+00) 

INTEGER (NYF04, NYL04 ) ;  NYF04=6;  NYL04=6 
YFRAC (6) =  2.640E-01 

*  Equal  Grid  spacing  in  Subregion  5 

*  SUBGRD (Y , 7,7,  4.400E-02,  1.000E+00) 

INTEGER ( NYF05 , NYL05) ?  NYF05=7;  NYL05=7 
YFRAC ( 7 ) =  3 . 080E-01 

*  Equal  Grid  Spacing  in  Subregion  6 

*  SUBGRD (Y , 8,9,  8.800E-02,  1.000E+00) 

INTEGER (NYF06, NYL06 ) ;  NYF06=8;  NYL06=9 
YFRAC (8 )  =  3.520E— 01; YFRAC (9 )  =  3.960E-01 

*  Equal  Grid  spacing  in  Subregion  7 

*  SUBGRD (Y , 10 , 11,  8 . 800E-02 ,  1.000E+00) 

INTEGER (NYFO 7 , NYL07 ) ;  NYF07=10;  NYL07=11 
YFRAC (10)=  4 . 400E-01 ; YFRAC ( 11) =  4.840E-01 

*  Equal  Grid  Spacing  in  Subregion  8 

*  SUBGRD (Y , 12,14,  1.320E-01,  1.000E+00) 

INTEGER (NYF08,  NYL08 )  ;  NYF08-^12  ;  NYL08  =  14 
YFRAC (12)=  5 . 280E-01 ; YFRAC ( 13 ) =  5.720E-01 
YFRAC (14 ) =  6 . 160E-01 

*  Equal  Grid  Spacing  in  Subregion  9 

*  SUBGRD (Y, 15,29,  6.600E-01,  1.000E+00) 

INTEGER (NYF09 , NYL09) ;  NYF09=15;  NYL09=29 
YFRAC ( 15) =  6 . 600E-01 ; YFRAC ( 16 ) =  7.040E-01 


69 


7 


YFRAC ( 17 )  = 
YFRAC ( 19 ) = 
YFRAC (21)= 
YFRAC (23)= 
YFRAC (25) = 
YFRAC (27)= 
YFRAC (29) = 


7 . 480E-01 ? 
8 . 360E-01 ; 
9. 240E-01? 
1. 012E+00; 
1. 100E+00; 
1. 188E+00; 
1 . 276E+00 


YFRAC (18) 
YFRAC (20) 
YFRAC (22 ) 
YFRAC ( 24 ) 
YFRAC (26) > 
YFRAC ( 2  8 ) : 


7 . 920E-01 
8.800E-01 
9. 680E-01 
1.056E+00 
1.144E+00 
1. 232E+00 


************************************************************ 

* 

*  GROUP  5.  Z -direction  grid  specification. 

* 


* 

******************************  4 *: **************************** 
* 

*  GROUP  6.  Body-fitted  coordinates  or  grid  distortion. 
BFC=T 

S£TPT(1, 1,1, 0,0,0) 

SETPT(1, 11, 1, 0, .22,0) 

SETPT(59, 11, 1, 1. 276, .22,0) 

SETPT(59, 1, 1, 1. 276, 0, 0) 

DOMAIN (1,59,1,1,1,1) 

SETLIN (XC, XF+LNI* (XL-XF) ) 

SETLIN ( YC, YF) 

SETLIN (ZC, ZF) 

DOMAIN (1,1,1,11,1,1) 

SETLIN (XC,XF) 

SETLIN (YC, YF+LNJ  * ( YL-YF) ) 

SETLIN (ZC, ZF) 

DOMAIN (1, 59, 11, 11, 1, 1) 

SETLIN (XC, XF+LNI* (XL-XF) ) 

SETLIN (YC, YF) 

SETLIN ( ZC, ZF) 

DOMAIN (59,59,1,11,1,1) 

SETLIN (XC, XF) 

SETLIN (YC, YF+LNJ* (YL-YF) ) 

SETLIN ( ZC , ZF) 

DOMAIN (1,59,1,11,1,1) 

MAGIC (T) 

SETPT (1,15,1,0, .308,0) 

SETPT(59, 15, 1, 1.276, .308,0) 

SETPT (21, 15, 1, .415, .308,0) 

DOMAIN (1, 1, 11, 15, 1, 1) 

SETLIN (XC,XF) 

SETLIN (YC, YF+LNJ* (YL-YF) ) 

SETLIN (ZC, ZF) 

DOMAIN (59, 59, 11, 15, 1,1) 

SETLIN (XC, XF) 

SETLIN (YC, YF+LNJ* (YL-YF) ) 

SETLIN (ZC , ZF) 

DOMAIN (1,21, 15, 15, 1, 1) 
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SETLIN (XC, XF+LNI* (XL-XF) ) 

SETLIN (YC, YF) 

SETLIN ( ZC , ZF) 

DOMAIN (21,59,15,15,1,1) 

SETLIN (XC,XF+LNI* (XL-XF) ) 

SETLIN (YC,  YF) 

SETLIN (ZC, 2F) 

DOMAIN (21 ,2 1,11 ,15,1,1) 

SETLIN (XC,XF+LNJ* (XL-XF) ) 

SETLIN (YC, YF+LNJ  * (YL-YF) ) 

SETLIN (ZC,2F) 

DOMAIN (1,21,11,15,1, 1) 

MAGIC (T) 

DOMAIN (2 1,59, 11, 15, 1,1) 

MAGIC (T) 

SETPT(1, 59, 1,0, 1.276,0) 

SETPT( 59, 59, 1,1. 276, 1.276,0) 

DOMAIN (1,1,15,59,1,1) 

SETLIN (XC,XF) 

SETLIN (YC, YF+LNJ* (YL-YF) ) 

SETLIN (ZC, ZF) 

DOMAIN (1,59,59,59,1,1) 

SETLIN (XC,XF+LNI* (XL-XF)  ) 

SETLIN (YC, YF) 

SETLIN ( ZC, ZF) 

DOMAIN (59, 59, 15,59, 1,1) 

SETLIN (XC, XF) 

SETLIN (YC, YF+LNJ* (YL-YF)  ) 

SETLIN (ZC, ZF) 

DOMAIN (1,59,15,59,1,1) 

MAGIC (T) 

* 

* 

************************************************************ 

* 

*  GROUP  7.  Variables  stored,  solved  &  named. 

* 

*  Solve  for  the  PRESSURE 

*  (Slab-by-slab  Method)  *  (Arithmetic  Averaging) 

SOLVE (PI) 

*  Solve  for  the  X-DIRECTION  VELOCITY  COMPONENT 

*  (Slab-by-Slab  Method)  *  (Arithmetic  Averaging) 

SOLVE (U1 ) 

*  solve  for  the  Y-DIRECTION  VELOCITY  COMPONENT 

*  (Slab-by-Slab  Method)  *  (Arithmetic  Averaging) 

SOLVE (VI) 

* 
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a*************************'*********************************/ 

* 

*  GROUP  8.  Tetms  (in  differential  equations)  &  devices. 

* 

* 

************************************************************ 

* 

*  GROUP  9.  Properties  of  the  medium  (or  media). 

* 

*  Set  First-Phase  Density  Value 
RH01=  1 . 100E+00 

*  Set  Laminar  Kinematic  Viscosity  Value 
ENUL=  1 . 000E-04 

*  Select  K-E  Turbulen  i-'odel 

*  ENUT  =  CMU  *  (Mixing  . jngth)  *  K**0.5 

*  ELI  =  (CD  *  K* *1.5) /E 
TURMOD ( KEMODL) 

* 

************************************************************ 

* 

*  GROUP  10.  lnter-phase-transfer  processes  and  properties 

* 

* 

*****************************  ****************************** 

* 

*  GROUP  ll.  Initialization  of  variable  or  porosity  fields 

* 

*  Initialize  the  X-DIRECTION  VELOCITY  COMPONENT 
FI INIT (U1 ) *  2 . 4  00E+00 

*  Initialize  the  TURBULENT  KINETIC  ENERGY 
FIINIT(KE) *  ° . 060E— 01 

*  Initialize  the  KINETIC-ENERGY  DISSIPATION  RATE 
FI INIT (EP)*  1 . 230E-01 

*  Obstructed  Region,  Number  1 

CONPOR (0.0, west , -NXF0 3 , -NXL03 , -NYF04 , -NYL05 ,1,1) 

*  Obstructed  Region,  Number  2 

CONPOR ( 0 . 0 , cell , -NXF0 1 , -NXL01 , -NYF01 , -NYL03 ,1,1) 

* 

*********************44*******4*4*4r***6*‘>  ******************* 

* 

*  GROUP  12.  Convection  and  diffusion  adjustments. 


**  * 


********************************  ************************ 


*  GROUP  13.  Boundary  conditions  and  special  sources. 

* 

*  INLET  Boundary  Condition,  Named  A 
PATCH (A, WEST, NXF01 , WXF01 , NYF04 , NYL04 , 1, 1, 1, 1) 

COVAL(A,  PI  ,  FIXrlU,  RHOl*  2.3',0E+0  >) 

COVAL ( A , U 1 , 0  .  .  YNS,  2. 360E+' 


COVAL ( A ,  VI , ONLYMS ,  O.OOOE+OO) 

COVAL ( A , KE , ONLYMS ,  1.880E-01) 

COVAL ( A ,  E  P , ONLYMS ,  1.070E-01) 

*  INLET  Boundary  Condition,  Named  C 
PATCH(C,WEST,NXF01,NXF01,NYF05,NYL05, 1, 1, 1, 1) 

COVAL (C, PI , FIXFLU, RHOl*  2.470E+00) 

COVAL (C,U1, ONLYMS,  2.470E+00) 

COVAL (C, VI, ONLYMS,  O.OOOE+OO) 

COVAL(C,KE, ONLYMS,  2.060E-01) 

COVAL (C,EP, ONLYMS,  1.230E-01) 

*  INLET  Boundary  Condition,  Named  D 
PATCH ( D , WEST , NXFO 1 , NXFO 1 , NYF06 , NYLO  ? , 1 , 1 , 1 , 1 ) 

COVAL (D, PI, FIXFLU, RHOl*  2 . 560E+00) 

COVAL (D,U1, ONLYMS,  2.560E+00) 

COVAL(D, VI, ONLYMS,  O.OOOE+OO) 

COVAL(D,KE, ONLYMS,  2.210E-01) 

COVAL(D,EP, ONLYMS,  1.370E-01) 

*  INLET  Boundary  Condition,  Named  E 
PATCH (E , WEST , NXF01 , NXF01 , NYF08 , NYLO 9 , 1 , 1 , 1 , 1) 

COVAL (E, PI, FIXFLU, RHOl*  2 . 560E+00) 

COVAL (E,U1, ONLYMS,  2.560E+00) 

COVAL (E, VI, ONLYMS,  O.OOOE+OO) 

CO VAL(E,KE, ONLYMS,  2.210E-01) 

COVAL(E,EP, ONLYMS,  1.370E-01) 

*  OUTLET  Boundary  Condition,  Named  OUTLET 
PATCH (OUTLET , EAST , NXL04 , NXL04 , NYF01 , NYL09 , 1 , 1 , 1 , 1 ) 
COVAL (OUTLET , PI , FIXP,  O.OOOE+OO) 

COVAL (OUTLET , KE , ONLYMS ,  0 . 000E+00) 

COVAL(OUTLET,EP, ONLYMS,  O.OOOE+OO) 

*  WALL  Boundary  Condition,  Named  TOPWALL 

PATCH (TOPWALL, NWALL, NXF01 , NXL04 , NYL09 , NYL09 , 1 , 1 , 1 , 1 ) 
COVAL (TOPWALL, U1,GRND2,  O.OOOE+OO) 

COVAL ( TOPWALL , KE , GRND2 , GRND2 ) 

COVAL (TOPWALL, EP, GRND2 , GRND2 ) 

*  WALL  Boundary  Condition,  Named  STEPWAL1 

PATCH (STEPWAL1 , NWALL, NXF01 , NXL01 , NYL03 , NYL03 , 1 , 1 , 1 , 1) 
COVAL (STEPWAL1 , U1 , GRND2 ,  0 . OOOE+OO) 

COVAL (STEPWAL1 , KE , GRND2 , GRND2 ) 

COVAL (STEPWAL1 , EP, GRND2 , GRND2 ) 

*  WALL  Boundary  Condition,  Named  STEPWAL2 

PATCH (STEPWAL2 , EWALL, NXL01 , NXL01 , NYF01 , NYLO 3 , 1 , 1 , 1 , 1) 
COVAL ( STEPWAL2 , V 1 , GRND2 ,  O.OOOE+OO) 

COVAL (STEPWAL2 , KE, GRND2 , GRND2 ) 

COVAL (STEPWAL2 , EP, GRND2 , GRND2 ) 

*  WALL  Boundary  Condition,  Named  BOTWALL 

PATCH ( BOTWALL, SWALL, NXFC2 , NXL04 , NYF01 , NYF01 , 1,1,1, 1) 
COVAL ( BOTWALL, U1 , GRND2 ,  0 . OOOE+OO ) 

COVAL ( BOTWALL, KE , GRND2 , GRND2 ) 

COVAL ( BOTWALL, EP, GRND2 , GRND2 ) 

* 
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************************************************************ 


★ 


*  GROUP  14.  Downstream  pressure  for  PARAB>= .  TRUE . . 

* 


* 

************************************************************ 


* 


*  GROUP  15.  Termination  of  sweeps. 

* 


*  Number  of  Iterative  Sweeps 
LSWEEP=600 

*  Reference  Residual  for  the 
RESREF (PI ) =  1 . 000E-04 

*  Reference  Residual  for  the 
COMPONENT 

RESREF (U1 ) =  1 . 000E-04 

*  Reference  Residual  for  the 
COMPONENT 

RESREF ( VI ) =  1.000E-04 

*  Reference  Residual  for  the 
RESREF (KE) =  1.000E-04 

*  Reference  Residual  for  the 
RATE 

RESREF (EP)=  1 . OOOE-04 

* 


(Outer  Iterations) 

PRESSURE 

X-DIRECTION  VELOCITY 

Y-DIRECTION  VELOCITY 

TURBULENT  KINETIC  ENERGY 
KINETIC-ENERGY  DISSIPATION 


************************************************************ 

* 

*  GROUP  16.  Termination  of  iterations. 

* 


* 

★ 

*  GROUP  17.  Under-relaxation  devices. 

* 

*  Automatic  False-Time-Step  Relaxation  Applied  to  Ui 
REAL(SCALEL, SCALEU) ;SCALEL=  4 . 400E-02 ;SCALEU=  2.560E+00 
RELAX (Ul, FALSDT,  1 . 500E+00*SCALEL/SCALEU) 

*  Automatic  False-Time-Step  Relaxation  Applied  to  Vi 
RELAX (VI # FALSDT,  1 . 500E+00*SCALEL/SCALEU) 

*  Automatic  False-Time-Step  Relaxation  Applied  to  KE 
RELAX (KE , FALSDT,  1 . 600E+00*SCALEL/SCALEU) 

*  Automatic  False-Time-Step  Relaxation  Applied  to  EP 
RELAX (EP, FALSDT,  1 . 600E+00*SCALEL/SCALEU) 

* 
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************************************************************ 

* 

*  GROUP  18.  Limits  on  variables  or  increments  to  them. 

* 

* 

************************************************************ 

* 

*  GROUP  19.  Data  communicated  by  satellite  to  GROUND. 

* 

* 

************************************************************ 

* 

*  GROUP  20.  Preliminary  print-out. 

* 

*  Deactivate  Print-Out  of  Satellite  Data. 

ECHO=F 

* 

★★★★a******************************************************* 

★ 

*  GROUP  21.  Print-out  of  variables. 

* 

*  Printout  for  the  PRESSURE 
OUTPUT (P1,Y,N,N,Y,Y,Y) 

*  Printout  for  the  X-DIRECTION  VELOCITY  COMPONENT 
OUTPUT (U1/Y,N,N,Y<Y,Y) 

*  Printout  for  the  Y-DIRECTION  VELOCITY  COMPONENT 
OUTPUT (VI, Y,N,N,Y,N,N) 

*  Printout  for  the  TURBULENT  KINETIC  ENERGY 
OUTPUT (KE,Y,N,N,Y,Y,Y) 

*  Printout  for  the  KINETIC-ENERGY  DISSIPATION  RATE 
OUTPUT (EP,Y,N,N,Y,Y,Y) 

* 

************************************************************ 

* 


*  GROUP  22.  Spot-value  print-out. 

* 

*  X-Direction  Index  of  Spot-Value 
IXM0N=6 

*  Y-Direction  Index  of  Spot-Value 
IYM0N=3 

* 

************************************************************ 

* 

*  GROUP  23.  Field  print-out  and  plot  control. 

* 

*  Print  TABLES  AND  PLOTS  of  Spot-Values  and  Residuals 
ITABL-3 

* 
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************************************************************ 

* 

*  GROUP  24.  Preparations  for  continuation  runs. 

* 

*  Save  the  Final  Flow  Field  in  a  RE-START  File 
SAVE=T ;  NSAVE=CHAM 

* 

************************************************************ 

STOP 
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